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The main purpose of the lecture was to lead students and young post-docs to 
the frontier of the hydrodynamic description of relativistic heavy-ion collisions 
(H.I.C.) in order for them to understand talks and posters presented in the 
Quark Matter 2008 (QM08) conference in Jaipur, India [1]. So the most recent 
studies were not addressed in this lecture as they would be presented during 
^ ' the QM08 conference itself. Also, we try to give a very pedagogical lecture 

^ . here. For the readers who may want to study relativistic hydrodynamics and 

its application to H.I.C. as an advanced course, we strongly recommend them 
to consult the references. 

^ ' This lecture note is divided into three parts. In the first part we give a 

. brief introduction to relativistic hydrodynamics in the context of H.I.C. In 

the second part we present the formalism and some fundamental aspects of 
relativistic ideal and viscous hydrodynamics. In the third part, we start with 
some basic checks of the fundamental observables followed by discussion of 
collective flow, in particular elliptic flow, which is one of the most exciting 
' phenomenon in H.I.C. at relativistic energies. Next we discuss how to for- 

, mulate the hydrodynamic model to describe dynamics of H.I.C. Finally, we 

conclude the third part of the lecture note by showing some results from ideal 
hydrodynamics calculations and by comparing them with the experimental 
^ ' data. 

^ , We use the natural units c = fi ~ ks — ^ and the Minkowski metric 

5^" = diag(l, —1, —1, —1) throughout the lecture note. 



2 Introduction to hydrodynamics in relativistic 
heavy-ion collisions 

The excitement raised by the announcement of the discovery of the "perfect" 
liquid at Relativistic Heavy Ion Collider (RHIC) in Brookhaven National Lab- 
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oratory (BNL) [2] is based on an agreement between predictions from ideal 
hydrodynamic models with the experimental data. While this agreement was 
certainly a large boost for various groups around the world doing research 
in hydrodynamics (and even in string theory!), there are also other reasons 
why the usage of hydrodynamics is strongly needed in H.I.C. Needless to 
say, the main goals of the physics of H.I.C. are to discover the deconfined 
nuclear matter under equilibrium, namely the Quark Gluon Plasma (QGP), 
and to understand its properties such as equation of state (EoS), tempera- 
ture and order of phase transition, transport coefficients and so on. However, 
the system produced in H.I.C. dynamically evolves within time duration of 
the order of 10-100 fm/c. Hence one has to describe space-time evolution of 
thermodynamic variables to fill a large gap between the static aspects of QGP 
properties and dynamical aspects of H.I.C. It is the hydrodynamics that plays 
an important role in connecting them. Various stages of H.I.C. are depicted 
in Fig. 1. Two energetic nuclei are coming along light-cone and collide with 
each other to create a multi-parton system. Through secondary collisions the 
system may reach thermal equilibrium and the QGP can be formed. This is 
a transient state: After further expansion and cooling the system hadronizes 
again. Eventually, expansion leads to a free-streaming stage through freezeout 
and particle spectra at this moment are seen by the detector. Hydrodynam- 
ics is applied to matter under local equilibrium in the intermediate stage. Of 
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Fig. 1. A schematic view of dynamics of a heavy ion collision along the collision 
axis. 

course, it is non-trivial a priori whether one can always apply hydrodynamics 
to the dynamics of H.I.C. Nevertheless it is not a bad idea to dare to apply it 
since we eager to understand the matter under equilibrium in terms of H.I.C. 

There is also another good reason to apply hydrodynamics to H.I.C. A lot 
of experimental data have been published so far at various collision energies. 
Ideally, one may want to describe these data from the first principle, i.e., 
quantum chromo dynamics (QCD). The QCD Lagrangean density reads 
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(1) 



where '04 is a quark field, D'' is a covariant derivative, m is a quark mass 
and Fj^" is a field strength of gluons. However, in spite of its simple-looking 
Lagrangean, it is very difficult to make any predictions directly from QCD 
in H.I.C. due to its complexity which mainly arises from non-linearity of in- 
teractions of gluons, strong coupling, dynamical many body system and color 
confinement. One promising strategy to connect the first principle with phe- 
nomena is to introduce hydrodynamics as a phenomcnological theory. We call 
this strategy a bottom-up approach to H.I.C. An input to this phenomcno- 
logical theory comprises the equation of state, 



which expresses the pressure P as a function of energy density e and baryon 
density n. This can be obtained by exploiting lattice numerical simulations of 
QCD'' . In the case of viscous hydrodynamics we need additionally the trans- 
port coefficients such as shear viscosity ry, bulk viscosity heat conductivity 
A, ctc.^ 

Once hydrodynamics turns out to work quite well in description of dynam- 
ics, one can utilize its outputs such as local temperature or energy density for 
other observables. In the current formalism of jet quenching, one needs an 
information of parton density or energy density along a trajectory of an en- 
ergetic parton [3, 4]. If one assumes J /t}) melts away above some temperature 
[5], one needs local temperature at the position of J /ip. In the case of electro- 
magnetic probes [6, 7], one convolutes emission rate (the number of produced 
particles per unit space-time volume at temperature T) of thermal photons 
and dileptons over the space-time volume under equilibrium. Hydrodynamics 
provides us with the information of the bulk matter. Therefore we can say 
that, in the context of H.I.C, hydrodynamics is the heart of the dynami- 
cal modeling: It not only describes expansion and collective flow of matter 
but also provides important informations in the intermediate stage for other 
phenomena. 

3 Formalism of the relativistic ideal/ viscous 
hydro dynamics 

The second part of the lecture note is more formal with many equations, 
but we try as much as possible to provide the intuitive picture behind the 

^ Prom lattice calculations, pressure as a function of temperature rather than energy 
density is obtained. Note also that, due to sign problem, thermodynamic variables 
are available only near the region of vanishing chemical potential. 

In principle, the information about these quantities can be obtained also from 
the lattice QCD simulations although it is much harder than the EoS. 



P = P(e,n), 



(2) 
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equations. The following references might be very helpful to complement this 
section [8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24]. 

3.1 The basic equations 

The basic hydrodynamical equations are energy-momentum conservation 

9pT'"'-0, (3) 
where T'^'^ is the energy-momentum tensor, and the current conservation 

d^Nt = i), (4) 

where Nj^ is the i-th conserved current. In H.I.C., there are some conserved 
charges such as baryon number, strangeness, electric charges and so on. We 
mainly assume the not baryon current N^g as an example of iV/* in the fol- 
lowing. In the first step we decompose the energy-momentum tensor and the 
conserved current as follows: 

T"'' = euV - PAf"" + WV + Wu^" + tt"^ , (5) 
TVf = riiu'' + y.'* . (6) 

All the terms in the above expansion will be discussed one by one later. Now 
we indicate that is the time-like, normalized four-vector 

u^u^ = 1 , (7) 

while the tensor A>^" is defined in the following way, 

Zi"" = - , (8) 

where g^^ is the Minkowski metric. We refer to and A^^ as the "projec- 
tion" vector and tensor operators, respectively. In particular, is the local 
flow four-velocity, but a more precise meaning will be given later, is per- 
pendicular to A^'^ ^ as can easily be seen from the definition of A^" given in 
Eq. (8) and from the fact that is normalized, 

= u^((7''" - u'^u") = - 1 • = . (9) 

Next we define the local rest frame (LRF) as the frame in which has only 
the time-like component non- vanishing and in which A^^ has only the space- 
like components non- vanishing, i.e., 

<RF = (1.0,0,0), (10) 
zij^^p =diag(0, -1,-1,-1). (11) 

As is easily understood from the above equations, one can say that u^{A'^'^) 
picks up the time-(space-)like component(s) when acting on some Lorentz 
vector /tensor. 

We now disc;uss the physical meaning of each term in the expansion of the 
energy- momentum tensor (5) and the conserved current (6). 
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Decomposition of T^^^ 

The new quantities which appear on the RHS in the decomposition (5) are 
defined in the following way: 

e = M^r'"^u^ (energy density) , (12) 

P = Ps + n = -^Af.uT^'" (hydrostatic + bulk pressure) , (13) 

W = A''^T°'''^U(j (energy (or heat) current) , (14) 
tt'"' = {T'"') (shear stress tensor) . (15) 

Each term corresponds to projection of the energy momentum tensor by one 
or two projection opcrator(s), and A^^"^. The first two equalities imply that 
the energy density e can be obtained from the time-like components of the 
energy- momentum tensor, while the pressure P is obtained from the space- 
like components. Contracting the energy- momentum tensor simultaneously 
with u'^ and Z\^'' gives the energy (heat) current W^. Finally, the angular 
brackets in the definition of the shear stress tensor tt^'' stand for the following 
operation, 



\{A^^A^p + A^^A\)-\a^^A^P 



A"!^ . (16) 



This means that {A^^") is a symmetric and traceless tensor which is transverse 
to w'' and . More concretely, one can first decompose the energy momentum 
tensor by two projection tensors symmetrically 

j^i^'^ = ^{A^^T^'f^Af^'' + A'^^T'^^A^") (17) 

and then decompose it once more into the shear stress tensor (traceless) and 
the pressure (non-traceless) 

(18) 



Decomposition of JV** 

In the decomposition (6) we have introduced the following quantities, 

Tii = u^N^ (charge density) , (19) 
= Zi'^^TV^ (charge current) . (20) 

The physical meaning of rij and Vj^ is self-evident from the properties of 
projection operators. 
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QUESTION 1: Count the number of unknowns in the above decompositions 
and confirm that it is IQiTi^") + Ak{N^) . Here k is the number of independent 
currents^. 



The various terms appearing in the decompositions (5) and (6) can be grouped 
into two distinctive parts, which we call ideal and dissipative part. In partic- 
ular, for the energy momentum tensor we have, 

T^"^ = Ti;" + ST'''' , (21) 
To"" = eu^u" - PsAf"" , (22) 
ST"" = -HA'"' + WV + Wu" + -k"" , (23) 

while for one charge current we have, 

TV' = + ^ATM , (24) 
= nu'' , (25) 
SN'' = V . (26) 

In the above relations Tl^''{N^) denote the ideal part, while the ST''''{5Nt') 
denote the dissipative part of the T'^''{N'^). 



3.2 The meaning of w** 

As we have already mentioned in Sec. 3.1, is the four-velocity of "flow". 
Now we would like to clarify what kind of flow we have in mind in this de- 
scription. In literature two definitions of flow can be found, 



1. flow of energy (Landau) [10]: 



T''„ul 1 



2. flow of conserved charge (Eckart) [8]: 



, (27) 
e 



(28) 



(see Fig. 2)^. In the first definition, also appears in the RHS of Eq. (27) 



If you consider as independent variables, you need additional constraint for 
them since these are redundant ones. If you also consider Pg as an independent 
variable, you need the equation of state Ps = Ps{e,n). 

Other definitions can be made. The situation here is quite similar to the gauge fix- 
ing condition iu gauge theories to eliminate the redundant variables. An essential 
point is to choose some "gauge" for later convenience. 
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Fig. 2. A sketcii of Landau and Ecltart definitions of flow. Two boxes are fluid 
elements. There is a "leak" current or according to the definition of fiow. 



So it should be understood as an equation with respect to m^. One may 
solve an eigenvalue problem for a given energy-momentum tensor . is 
a normalized time-like eigenvector and the corresponding positive eigenvalue 
is energy density e. If the dissipative currents are small enough, one can show 
the following relation between these two definitions of flow 



(29) 



Equation (29) can be shown by assuming that both two definitions of flow can 
be connected by infinitesimal proper Lorentz transformation 



(J^ 



(30) 
(31) 



where e^^ infinitesimal anti-symmetric tensor and neglecting the higher orders 
of dissipative currents. Obviously, = {V^ — 0) in the Landau (Eckart) 
frame. In the case of vanishing dissipative currents, both definitions represent 
a common flow. In other words, flow is uniquely determined in the case of 
ideal hydrodynamics. We should emphasize that Landau definition is more 
relevant in the context of H.I.C. at ultrarelativistic energies since we expect 
a small baryon number is deposited near the midrapidity region. 



3.3 Entropy 

We start this subsection by briefly discussing the entropy conservation in 
"ideal hydrodynamics" . By "ideal hydrodynamics" we mean the case when en- 
tropy is not produced during the evolution^. Neglecting the dissipative parts, 
the energy-momentum conservation (3) and the current conservation (4) re- 
duce to 

Note that, if discontinuities exist in the solution, entropy is produced even in 
ideal hydrodynamics. 
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W = 0, (32) 
W = 0, (33) 

where T//'^ and are the ideal parts introduced in Eqs. (22) and (25). 
Equations (32) and (33) are the basic equations of "ideal hydrodynamics". 
By contracting Eq. (32) with it follows, 

= u^d^To^'' 

= T(M^a^s + sdi^u'') + niui'd^n + n^^u^) . (34) 

We have introduced here the temperature T , entropy density ,s and chemical 
potential /i through the first law of thermodynamics de — Tds + ^dn. Here it 
is assumed that tliermalization is maintained locally. The second term on the 
RHS in Eq. (34) vanishes due to Eq. (33). If we now introduce the entropy 
current as 

5" = su" , (35) 

it follows from Eq. (34) that 

df.S" = d^isu'') = ^i^'d^.s + sd^uf" = . (36) 
Hence the entropy is conserved in ideal hydrodynamics. 

QUESTION 2: Go through all steps in the above derivations. 

Now we go back to viscous hydrodynamics. Hereafter we consider only the 
Landau frame and omit the subscript L. For simplicity, we further assume that 
there is no charge in the system although in the realistic case a small amount 
of charge might exist in the system. What we are constructing here is the 
so-called first order theory of viscous hydrodynamics. The main assumption 
is that the non-equilibrium entropy current vector S'^ has linear dissipative 
term(s) constructed from V^, U and n'^" and can be written as 

S-^ = su" + aV . (37) 

The first term on the RHS is the ideal part and the second term is the cor- 
rection due to the dissipative part. It is impossible to construct a term which 
would form a Lorentz vector from w^" on the RHS in the above equation be- 
cause tt'"' is perpendicular to 77/* by definition*. Since we have also assumed 
that there is no charge in the system, i.e., = 0, it follows that aV^ van- 
ishes. 



Also remember = in the Landau definition. One may think that 77m'' 
is a possible candidate in the entropy current S'^. However, the second law of 
thermodynamics is not ensured in this case. See also discussion in Ref. [21]. 
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We now calculate the product of the temperature T and the divergence of 
the entropy current (37). It follows, 

= TT^, {V^u") - nd^u^" . (38) 

where V = A^^^d^. In transferring from the second to third line in the above 
calculation we have used the energy-momentum conservation, d^T'^'' = 0. It 
is very important to note that due to the assumption that there is no charge 
in the system we could neglect the dissipative part of entropy current (37), 
but the dissipative part of energy-momentum tensor (23) does not vanish. The 
non-vanishing dissipative part of energy-momentum tensor gives a contribu- 
tion which yields a difference between the equations characterizing the first 
order theory of viscous hydrodynamics and the equations of ideal hydrody- 
namics derived before. 

QUESTION 3: Check the above calculation. 

In order to solve the hydrodynamic equations we must first define the dissi- 
pative current. We introduce the following two phenomenological definitions, 
so-called constitutive equations, for the shear stress tensor 'k'^" and the bulk 
pressure iT, 

it'"' = 2r] {V^u") , (39) 
n = -Cd^ui" = -CV^w" . (40) 

In Table 1 we outline the new variables and terminology used in the above 
equations. Notice that, within our approximation A'^'^ = 0, there is no vector 
component of thermodynamics force. 



Table 1. New variables and terminology. 



Thermodynamic force 


Transport coefficient 


Current 


X"" = (Vm") 
tensor 


V 

shear viscosity 




X = -a^w'^ 
scalar 


c 

bulk viscosity 


n 



After inserting the definitions (39) and (40) in the last line of (38), we arrive 
at, for positive transport coefficients. 
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^"^^^ - 277 C 

= 27?(W)' + C(-a^M''f >0. (41) 

This ensures the second law of thermodynamics 

9^5'^ > . (42) 

In the case of viscous hydrodynamics, entropy is not decreasing. 

3.4 The equations of motion 

In order to derive the equations of motion, we use again energy-momentum 
conservation (3). After contracting Eq. (3) with Ui, we have 

u^d^.T'"' = , (43) 

from which one can obtain the first equation of motion, 

e = -{e + Ps+n)0 + ni,, {V^u") . (44) 

On the other hand, after contracting Eq. (3) with it follows, 

^^a9^T«'5 = o, (45) 

from which one can obtain the second equation of motion, 

(e + + n)uf = Vt\Ps + 77) - A^"^Vf^^lap + 7r^"M„ . (46) 

This is exactly the relativistic extension of the Navier-Stokes equation. In 
writing the above equations we have introduced, 

= d^u^ expansion scalar (divergence of flow), (47) 
"dot" = D = u^d'^ substantial time derivative . (48) 



QUESTION 4: Starting from the energy-momentum conservation (3) derive 
equations (44) and (46). 

To get some intuitive interpretation of the first equation of motion, we insert 
expressions (39) and (40) for the shear stress tensor and bulk pressure into 
Eq. (44), 

e = -ee-Ps9+ti- + ^-^ 

C 2r] 

= -ee-P,e + C{-Of + 2r]{V''u''f . (49) 
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The above equation determines the time evolution of energy density e in the 

CO- moving system. The first term on the RHS describes dilution/compression 
of energy density due to the change of volume, because can be expressed in 
terms of volume of a fluid element V as 

0^^. (50) 

In ideal hydrodynamics, this relation holds exactly. If the system expands 
{0 > 0), the energy density is diluted. So the effect of expansion appears as 
negative source term —e0 in Eq. (49). If we move along with a fluid element, 
the internal energy in the fluid element is not conserved due to the work 
done by pressure, which is described by the second term on the RHS in (49). 
Finally, the last two positive definite terms in (49) represent the production 
of entropy which heats up the system. 

Now we comment on the second equation of motion (46). But before doing 
that, we recall the non-relativistic Navier-Stokes equation, 

Dv = --VPs + -V^v. (51) 
P P 

Here p is the mass density, r] is shear viscosity and D = ^ + v • V is the non- 
rclativistic version of substantial time derivative. The above version of the 
non-relativistic Navier-Stokes equation applies to the case of incompressible 
fluids such that V • u = is valid. On the LHS we have the time derivative of 
velocity, which is nothing but acceleration. The first term on the RHS is the 
source of the fiow and it is solely due to the pressure gradient VP,, while the 
second term represents the diffusion of the flow. The flnal flow velocity comes 
from the interplay between these two terms: The first term generates the 
flow, while the second term dilutes it. The ratio rj/p is called kinetic viscosity 
and plays a role of diffusion constant in the Navier-Stokes equation (51). The 
diffusion term in Eq. (51) requires more detailed treatment. For an illustrative 
purpose, consider first the heat equation in (A^-l-l)-dimensional space-time 

i ' 

where T is temperature and constant k is heat conductivity in some unit. One 
can discretize the heat equation (52) in (2 + l)-dimensional space-time: 



rpn-\-l rpn I 



rpn I rpn I rpn I rpn 



rpn 



where i and j are indices of the site and n is the time stcip. The first term in the 
brackets in Eq. (53), Tj^-, indicates an average of temperature around the cell 
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under consideration. If temperature at (z, j) is smaller (larger) than the aver- 
aged one Ti j > Ti j (Ti j < Tij), the second term in Eq. (53) becomes positive 
(negative) and, consequently, temperature increases (decreases) in the next 
time step. Repeating this procedure, temperature becomes flat even if starting 
from a bumpy initial condition. Thus, generally speaking, the second deriva- 
tive with respect to coordinates describes averaging/smoothening/diffusion 
of given distributions and a coefiicient in front of it describes how quick the 
distribution diffuses. Now going back to the Navicr-Stokcs equation (51), it is 
obvious from the above discussion that the second term describes diS'usion of 
flow and that kinetic viscosity ij/p plays a role of a diffusion constant. The rel- 
ativistic version of Navier-Stokcs equation (46) has a similar form to Eq. (51) 
if one plugs in constitutive equations (39) and (40) and assumes the fluid is 
incompressible, ^ = 0. 

Bjorken's equation in the 1** order theory 

Now wc rewrite again the first equation of motion by making use of Bjorken's 
ansatz [16] 



where r = VW—z^ is the proper time. This is a boost invariant Bjorken's 
solution which is also called 1-dimensional Hubble flow since velocity in the 
z direction, v^, is proportional to z, which is an analogy to three dimensional 
Hubble flow of the universe. After inserting this solution into the constitutive 
equations (39) and (40) 




(54) 




(55) 



c 



= diag(l, 0,0,-1 



) 



(56) 



n 



(57) 



r 



we arrive at the following equation of motion 




(58) 



This equation determines the time evolution of energy density in the 1^* order 
theory in 1-dimensional expansion. 



QUESTION 5: Derive equation (58). 



On the RHS of (58) wc have three terms in the bracket. If wc neglect the last 
two terms this equation reduces to the famous Bjorken equation [16] which 
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states that in ideal hydrodynamics the energy density evolution is determined 
by the sum of energy density e and the hydrostatic pressure Pg , divided by the 
proper time r. The last two terms on the RHS in (58) represent the viscous 
correction to ideal hydrodynamics. The first one is the viscous correction 
originating from the shear viscosity in compressible fluids, while the second 
one comes from the bulk viscosity. We remark that both terms are proportional 
to 1/r which is due to the fact that the expansion scalar 6 in the Bjorken 
scaling solution can be written as 

e=-. (59) 

T 

Two transport coefficients in the viscous correction, rj/s and C,/s, turn out to 
be very important. They are the dimensionless quantities in natural units and 
reflect the intrinsic properties of the fluids^. 

Recently progress has been made in obtaining the values of the transport 
coefficients from microscopic theories. Here we summarize the most important 
results and conclusions, 

• rf/s = 1 1 An and C,/s = Q are obtained from Af = 4 SUSY Yang- Mills 
theory [25]. The latter one is automatically obtained from the conformal 
nature of the theory; 

• ri/s = 0(0.1 — 1) for giuonic matter is obtained from the lattice calculations 
of pure SU(3) gauge theory [26]; 

• Bulk viscosity has a prominent peak around Tc resulting from trace 
anomaly of QCD [27, 28] (see also a phenomenological approach in 
Ref. [29]). 

3.5 The 2"** order theory £ind its application to Bjorken's equation 

There is an important issue in the first order theory which is the violation 
of causality. We can trace back the origin of the violation of causality to 
our phenomenological definitions (39) and (40) for the shear stress tensor 
and the bulk pressure, respectively, and to the fact that the Navier-Stokes 
equation is a parabolic equation, namely, the time derivative is of first order, 
while the space derivative is of second order. The same arguments hold also 
for the violation of causality in relativistic hydrodynamics: It is known that, 
under linear perturbations on the moving background equilibrium state, the 
solutions are unstable and acausal [30] (for a more detailed discussion, see 
also a recent study in Ref. [31]). For an illustrative purpose, wc continue this 
discussion by analyzing the heat equation as an example of the parabolic 
equation in three dimensional space^°, 

^ We stress that in the context of H.I.C. the statement which is often used, "viscos- 
ity is smah" , is not precise. From the equations we have derived we see that the 
correct statement should be "viscous coefficients are small in comparison with 
entropy density". 
^° Again, we choose some unit to simplify the following equations. 
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dT A 



dt -"^^dxf 

The heat equation can be easily derived by combining the balance equation, 

dT dq' 
-dt=-d^^' 

together with the constitutive equation, 

dT 

(f = —K— — Fourier's law. (62) 

In the above equations T is the temperature, is the heat current and k is 
the heat conductivity. The above constitutive equation is purely phcnomeno- 
logical. Although we are here considering the non-relativistic equations, the 
general arguments and conclusions we write down are valid in the relativistic 
case as well. The heat equation (60) violates causality. It can be easily con- 
firmed that the Green's function of the heat equation (60), sometimes called 
heat kernel, is Gaussian 



G{x\t;x'o,to) = rrrexp 

[4TrK{t-to)]2 



(X* - 



(63) 



and the "long tail" of this Gauss function causes the violation of causality in 
the heat equation. This issue was heuristically resolved by Cattaneo in 1948 
[9] after an additional term on the LHS of the constitutive equation (62) was 
introduced "by hand" , 

da' dT 

In the modified constitutive equation we have a new constant which is often 
called the "relaxation time" . Correspondingly, the heat equation (60) is also 
modified, 

^^^ + -dt=^d^' c. = V'^/r.. (65) 

In the literature the above equation is known as a telegraph equation. While 
the original heat equation can be classified as a parabolic equation, the tele- 
graph equation belongs to the family of hyperbolic equations. Causality is not 
violated in Eq. (65) simply because we can now, by choosing the relaxation 
time Tr to be large, reduce the signal velocity Cg to values smaller than the 
speed of light c. 

In relativistic hydrodynamics the relaxation terms introduced above can 
be obtained by modifying the entropy current in the following way. 



S"* = su^ + 0{ST^'') + O {{5T^"'f) 



(66) 
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By including the quadratic dissipative terms we are starting to work within the 
framework of 2""^ order theory. The non-cquihbrium entropy current vector 
in the 2"'^ order theory has hnear + quadratic dissipative term(s) constructed 
from (y*, 77, tt^'^). Again, we demand the 2"*^ law of thermodynamics, d^S^ > 
0. Thus, quadratic dissipative terms modify the constitutive equations which 
now read, 

^^^^a^./3^^^ + tt''^ = 2rj {V^u") + ■■■ , (67) 
Tnfl + n = -Ca^u" + • • • . (68) 

When compared to the constitutive equations of the 1*^' order theory. (39) 
and (40), we see that in the 2"^^ order theory in each constitutive equation 
a relaxation term appears. Relaxation terms include and rn, which are 
the relaxation times. It is important to note that in the 2"^^ order theory the 
constitutive equations are no longer algebraic equations. As a consequence, 
dissipative currents become dynamical quantities like thermodynamical vari- 
ables. The constitutive equations with relaxation terms have been employed in 

recent viscous fluid simulations [32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44] 
11 

Finally, we outline the Bjorken's equation in the 2"*^ order theory, 
dn 4t] ttt^ nijT d 

"'^rf7 + "=37-^-— rf^;^' 
dn^^ C nrn ncT d T„ 

dr T T 2 dr Ql 

where 

TT = 7r°° - n'' . (72) 

It is easy to show that the above formulas reduce to the ones in the 1^* order 
theory if one takes t^- ^ and Tn 0. We remark here that, contrary to the 
1^* order theory, one needs to specify initial conditions for dissipative currents 
in the 2"*^ order theory. 



3.6 Summary 

Let us summarize the main points so far: 

• Hydrodynamics is a framework to describe the space-time evolution of 
matter under local thermal equilibrium; 

Some of the references here do not employ the same equations as mentioned 
here. There arc still some hot debates how to formulate the correct relativtstic 
equation of viscous fluids or which terms in the constitutive equations of the 2"'' 
order theory should be kept in the simulations. 
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• A naive extension of Navier-Stokes equation to its relativistic version, 
which is called the first order theory, has problems on instabilities and 
causality; 

• Relaxation terms are needed in the constitutive equations to resolve the 

above issues; 

• These terms naturally arise in the constitutive equations when the 2"^* 
order corrections of dissipative currents are considered in the entropy cur- 
rent. 



4 Applications 

In this section we apply the formalism of hydrodynamics to heavy-ion colli- 
sions. As already noted in Sec. 1, we do not argue recent analyses in terms 
of viscous hydrodynamics. We show only results from ideal hydrodynamic 
models. One can also consult recent other reviews of hydrodynamic models at 
RHIC which complement the present lecture note [45, 46, 47, 48, 49, 50, 51, 
52]. We start by discussing some basic tests of whether the system produced 
in H.I.C. can be described by thermodynamic quantities. Then we discuss 
collective flow and introduce ideal hydrodynamic models to describe the flow 
phenomena in H.I.C. Finally we show results from ideal hydrodynamic models 
and compare them with experimental data. 



4.1 Basic checks of observables at RHIC 



Rc!eent lattice QCD results show [53] the energy density as a function of the 
temperature suddenly increases by ~ 1 GeV/fm^ at the (pseudo-)critical tem- 
perature Tc ~ 190 MeV Above this temperature, the system is supposed to 
be in the deconfined QGP. The first check is whether the energy density pro- 
duced in H.I.C. is sufHcient to form a QGP. Phenomenologically, the energy 
density in H.I.C. can be estimated through Bjorken's formula [16] 

Here (rriT) is the mean transverse mass, y = ^ In ^^.p" is the rapidity, ^ is 

the number of particles per unit rapidity, r = \/t'^ — is the proper time 

and _R is an effective transverse radiiis. The energy density obtained above 
depends on the proper time since the system is supposed to expand in the 



Energy density increases with temperature rapidly but smoothly. So this is not a 
phase transition but a cross over in a thermodynamically strict sense. This is the 
reason why we call it pseudo-critical temperature here. 

This formula neglects the effect of -pdV work. If the system is kinetically equi- 
librated, the energy density should be larger than the value obtained by this 
formula [54, 55]. 
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longitudinal direction with the expansion scalar 6 = l/r. One can compare 

Bjorkcn's energy density to the energy density from lattice QCD simulations 
to see whether it is sufficient energy density to form a QGP. Figure 3 shows 




Fig. 3. eBjT versus the number of participants at three collision energies [56]. 



the PHENIX data on sbjT versus the number of participants at three collision 
energies [56]. If r is taken to be 1 fm/c, the Bjorken's energy densities at 
^/sNN = 130 and 200 GeV are well above the energy density at the transition 
region ^ 1 GeV/fm^. Therefore sufficient energy is deposited in the central 
rapidity region in H.I.C. at RHIC. However, attention should be paid to the 
interpretation. The above formula just counts the total measured energy di- 
vided by the volume of a cylinder. So the system is not necessary thermalized. 
In this sense, this is a necessary condition, not a sufficient condition, to form 
a QGP. 

The next basic check is whether the matter in H.I.C. reaches chemical 
equilibrium. Assuming thermal and chemical equilibrium, we can calculate 
the number density of a certain particle species 

^^^^^^^=2^^i exp[(£;.-M.)/T]±l- ^''^ 

Hi gives the number density of particle species z as a function of the temper- 
ature T and chemical potential jii. g is the degeneracy of the particle, p is 
the momentum and E is the energy. We further assume the measured parti- 
cle number is fixed at a certain temperature and chemical potential, which is 
called chemical freezeout. Then the average number of particles, (A^i), can be 
estimated by summing contribution from particles directly emitted from the 
system with volume V and contribution from resonance decays 



{Ni) = V 



(75) 



Here nf^ and nji are the number density of dirc;c;tly emitted particle i and 
resonance R, respectively. -Tr^i is the branching ratio of the resonance R 
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decaying into species i. When one looks at ratios of two particle numbers, 
the volume V is canceled out. Thus the particle ratios depend only on two 
parameters: the temperature T and the baryonic chemical potential fis- In 
Fig. 4, various combinations of the particle ratio observed at RHIC are fitted 
by two parameters [57]. We find a remarkably good fit to data with only these 
two parameter^^. At ^/sWN = 130 GeV, the temperature is fitted to be 176 



P''P 



A/A Zll fl/Qit-/7c*KVK*K7ii- p/;i- K 7h" (|)/h" Alh' =/hQ/jt"*10 ; 

^ -s^ 'St 



t:r STAR 
■ , PHENIX 
□ PHOBOS 
A BRAHMS 

^FI=^'30 GeV 

Model re-fit with ^li data 
T = 17eiVleV, u =41MeV 



t 



i. 



+ 



p/p K7K*K77i- p/7t<.""2/h"'50 



Model prediction for 

T = 177MeV, u =29MbV 



Braun-Munzinger et ai., PLB 518 (2001)41 



D. Magestro (updated Juiy 22, 2002) 



Fig. 4. Ratios of particle number at RHIC [57]. 



MeV which is close to the critical temperature from lattice QCD calculations. 
At the temperature, which we call chemical freezeout temperature T"^^, the 
system ceases to be in chemical equilibrium. So we expect the system reaches 
chemical equilibrium above T'^^. Again, one has to keep in mind that this is a 
necessary condition since even in e"'"e~ or pp collisions observed particle ratios 
are fitted reasonably well by using statistical models [58] . See also discussions 
in, e.g., Refs. [59, 60]. 

The last basic check is whether the matter reaches kinetic equilibrium. If 
we suppose a system in H.I.C. is in kinetic equilibrium, the pressure is built 
inside the system. The matter is surrounded by vacuum, so pressure gradient 
in outward directions generates collective flow and, in turn, the system ex- 
pands radially. The momentum distribution in kinetically equilibrated matter 
is isotropic. On the other hand, when the matter is moving at a finite ve- 
locity the momentum distribution is Lorentz boosted. This is illustrated in 
Fig. 5. If this kind of distortion in momentum distribution can be observed 
experimentally, one can obtain some information about kinetic equilibrium. 
Assuming each fluid element expands radially at radial flow velocity vt, the 
Pt spectra for pions and protons can be calculated by convoluting these dis- 



There are some additional parameters in the recent statistical models such as 
excluded volume correction, strangeness suppression factor and so on for a better 
description of the data. 
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(a) Isotropic case (b) Lorentz boosted in positive x 

direction 

Fig. 5. Fluid elements at rest and at a finite velocity in x direction. Momentum 
distribution in the latter case is distorted by Lorentz boost along x axis. 

torted momentum distributions over azimuthal direction (blast wave model 
[61, 62]). Here pr is the transverse momentum which is perpendicular to the 
collision axis. The green curves are results with T = 100 MeV and radial flow 
velocity vt = 0.5. On the other hand, the red curves are results with T = 160 
MeV and vanishing flow vt — 0. For light particles like pions, there is al- 
most no sensitivity to distinguish the two cases: Reduction of temperature is 
almost compensated by radial flow. However, in the case of heavier particles 
like protons, a clear difference can be seen between these two cases: There is a 
shoulder structure at low pT resulting from radial flow. This kind of spectral 
change is observed in H.I.C., as can be seen in Fig. 7. It shows the proton pT 
spectra for p+p (black) , d+ Au (pink) and Au+ Au (red) collisions obtained by 
STAR Collaboration [63]. For p+p and d+Au collisions the spectra have just a 
power-law shape. However, in Au-|-Au collisions, one sees a shoulder structure 
at low pt (< 1 GeV/c). This is consistent with a thermal plus boost picture 
and suggests that a large pressure could be built up in Au-|-Au collisions. 
One can fit the pr spectrum using a blast wave parametrization [61, 62] and 
obtains decoupling temperature T'^'^'^ and the mean collective flow velocity as 
a function of the centrality. Even for pp collisions these parameters are finite 
(see Fig. 8) [64], which indicates that a more sophisticated model would be 
needed to interpret the data. This kind of spectral change can also be seen in 
results from kinetic theories in which kinetic equilibrium is not fully achieved. 
Therefore it is indispensable to perform a systematic study based on a more 
sophisticated dynamical framework. 
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Fig. 6. pt spectra for pions and protons from a thermal plus boost picture. See 
text for details. 
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Fig. 7. Proton spectrum for pp (black), dAu (pink) and AuAu (red) collisions. 
Adopted from a presentation file by O. Barannikova at Quark Matter 2005, Bu- 
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Fig. 8. Fitted parameters in blast wave model calculations [64]. 



We have obtained the necessary conditions for studying the QGP: (1) 
The energy density can be well above the critical value which is predicted 
from lattice QCD simulations; (2) A chemical freezeout temperature extracted 
from particle ratios is close to pseudo-critical temperature which is again 
from lattice QCD simulations; (3) High pressure can be built up in H.I.C., 
which suggests the system reaches kinetic equilibrium. If one of them was not 
confirmed through these basic checks, one would not need to go next steps 
towards detailed studies of the QGP in H.I.C. 



4.2 Elliptic flov^r 

Before going to a detailed discussion on the hydrodynamic models, we dis- 
cuss collective flow, in particular, anisotropic transverse flow. Here "collective 
flow" is meant by the correlation between position of matter and direction of 
flow, which is not necessary to be hydrodynamically evolving matter. A good 
example has already appeared in the previous subsection. In the case of radial 
flow, velocity of expanding matter has a component parallel to the radial co- 
ordinate. Figure 9 shows a heavy ion collision in the reaction plane (left) and 
transverse plane (right). In such a collision a region of the locally equilibrated 
state can be created. In the transverse plane the overlap region has an almond 
like shape, so the region is anisotropic with respect to the azimuthal angle. 
The azimuthal momentum distribution can be expanded into a Fourier series 

15 

dN N 

— = — [1 + 2wi cos((/)) + 2v2 cos(2(/)) + ■■■], (76) 
am Ztt 



Jd(t)COs{n<j))^^ 
T' 



^» = . , , rf^ ^ = (cos(n<^)) . (77) 
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Here we suppose azimuthal angle is measured from reaction plane. Of course, in 
the experimental situations, the reaction plane is not known a priori. We will not 
go into details of how to find reaction plane experimentally. 
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(a) In the reaction plane (b) In the transverse plane 

Fig. 9. Illustration of a H.I.C. 

where (j) is the azimuthal angle of momentum and w„ are the Fourier coefii- 
cients of n-th harmonics [65] . Because of the symmetry around the y-axis the 
sine terms vanish. The first and the second harmonics, Vi and V2, are called 
directed and elliptic flow parameters, respectively. The first harmonic, wi, is 
illustrated in Fig. 10(a). Particles are emitted preferably, e.g., in the direc- 
tion of the large arrows in the reaction plane. Directed flow is significantly 
seen near the beam rapidity region, but vanishes near midrapidity due to 
symmetry of the collision geometry. The second harmonic, V2, is much more 




(a) First harmonic vi (b) Second harmonic V2 

Fig. 10. Anisotropic transverse flow 



relevant for studying matter around midrapidity in H.I.C. at relativistic ener- 
gies since spectators already fly away [66] , therefore a lot of efforts to measure 
V2 have been made at RHIC so far. One of the first observables was actually 
V2 measured by STAR Collaboration [67]. It is illustrated in Fig. 10(b). 

Elliptic flow is how the system responds to the initial spatial anisotropy 
[66, 68, 69, 70]. Suppose two extreme situations illustrated in Fig. 11. In the 
first case (see Fig. 11(a)) the mean free path among the produced particles 
is much larger than the typical size of the system. In this case the azimuthal 
distribution of particles does not depend on azimuthal angle on average due 
to the symmetry of the production process. The other extreme case is when 
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(a) Large mean free path 



(b) Small mean free path 



Fig. 11. Normalized azimuthal distribution dN/d<j) of a non-central H.I.C. 

the mean free path is very small compared to the typical system size (see 
Fig. 11(b)). In this case hydrodynamics can be applied to describe the space- 
time evolution of the system. The pressure gradient along the horizontal axis 
is much larger than along the vertical axis due to the geometry. So the collec- 
tive flow is enhanced along the horizontal axis rather than along the vertical 
axis and, in turn, the azimuthal distribution gets oscillated. The amplitude 
of this oscillation in the normalized azimuthal distribution describes exactly 
the elliptic flow parameter. In this way, the elliptic flow is generated by the 
spatial anisotropy of the almond shape due to multiple interactions among the 
produced particles. We have good opportunities to extract some information 
about the mean free path from the elliptic flow analysis. 

The eccentricity is a very important quantity to interpret elliptic flow 
phenomena. To quantify the initial almond shape, the following formula can 
be used 



e = 



(78) 



The brackets denote an average over the transverse plane with the number 
density of participants as a weighting function 



= / dxdy ■ 



f-part 



(79) 



This is sometimes called the standard eccentricity. If the system is elongated 
along the y-axis, the eccentricity is positive. In more realistic situations, the 
eccentricity fluctuates from event to event. This fluctuation of the initial ec- 
centricity [71, 72, 73, 74, 75, 76] is important to understand the elliptic flow 
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in the small system such as Cu+Cu collisions or peripheral Au+Au collisions. 
Figure 12 shows a sample event projected into the transverse plane from a 




Fig. 12. An example of participants (magenta) and spectators (yellow and orange) 
in a H.I.C. from a Monte Carlo Glauber model. Adopted from a presentation file by 
D. Hofman at Quark Matter 2006, Shanghai, China. 



Monte Carlo Glauber model. Participants are shown in magenta and spec- 
tators are in yellow and orange. In this case one could misidentify the tilted 
line as the reaction plane, while the true reaction plane is the horizontal axis 
(dashed line). The angle of the tilted plane with respect to the true reaction 
plane fluctuates event by event. Of course we cannot observe the true reac- 
tion plane from experimental data. On the other hand, an apparent reaction 
plane (tilted line in Fig. 12) is determined also by elliptic flow signal itself. 
Another definition, called the participant eccentricity, is much more relevant 
for quantifying the almond shape in the event by event basis 







^xJ f xy 


^part 
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= {y'}~ 


- {y? , 




= {xy} - 





(80) 

(81) 
(82) 
(83) 

Now the average {• ■ •} is taken over in a single event generated by a Monte 
Carlo Glauber Model. 

In the following, the important properties of elliptic flow are demonstrated 
through hydrodynamic/transport simulations of H.I.C. In hydrodynamic sim- 
ulations, the eccentricity is usually defined by weighting local energy density 
e(x, y) or local entropy density s{x^ y) in the transverse plane rather than the 
number density of participants ?^part(2^, 2/)- Figure 13 shows the eccentricity 
Ex and the momentum eccentricity 

Jdxdy{TS^-Tyy) 
" Jdxdy{T--+Tyy) ^ ' 
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as a function of the proper time from a hydrodynamic simulation assuming 
Bjorken scaling solution in the longitudinal direction and two different sets 
of the EoS [45]. Details of hydrodynamic models will be discussed later. The 




Fig. 13. The spatial eccentricity Sx and the momentum eccentricity Sp as a function 
of the proper time r in Au+ Au collisions at 6 = 7 fm [45] . Solid and dashed curves 
correspond to two different sets of the EoS. 
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Fig. 14. V2/ex as a function of impact parameter b [77]. 

spatial eccentricity Sx decreases as the system expands and the momentum 
anisotropy rapidly increases at the same time. So the spatial anisotropy turns 
into the momentum anisotropy. The momentum anisotropy Sp is created and 
saturates in the first several femtometers, so the observed V2 is expected to 
be sensitive to the initial stage of the collision. Figure 14 shows the impact 
parameter dependence of the ratio of output (^2) to input (ex) [77] which can 
be understood as a response of the system. Ideal hydrodynamics predicts that 
V2 is roughly proportional to the eccentricity 




V2 « 0.2e 



(85) 



26 



Tetsufumi Hirano, Naomi van der Kolk, and Ante Bilandzic 



Figure 15 shows a result from a kinetic approach based on the Boltzmann 
equation for gluons undergoing elastic scattering only [78] Starting with 



1 mb 

free streaming 



Fig. 15. 112 as a function of proper time from Boltzmann calculations for different 
gluon cross sections [78]. Curves are guide to eyes. 



a uniform distribution in an almond shape in coordinate space and thermal 
distribution in momentum space, the multi-giuon system expands according 
to the Boltzmann equation with various transport cross sections^^. From this 
figure we can understand several important features of the elliptic flow: 

1. V2 is not generated in the free streaming case, so elliptic flow is generated 
indeed through secondary collisions; 

2. Elliptic flow is generated in the early stage of the collision and saturates 
after the first 2 to 3 fm/c; 

3. The saturated value of V2 is sensitive to the cross section among the par- 
ticles ^ ^ 

CTtr oc - cx - , (86) 
A 77 

where A is the mean free path and rj is the shear viscosity calculated in 
the kinetic theory of gases; 

4. In the limit of large transport cross sections (strongly interacting limit), 
the system is expected to reach the ideal hydrodynamic result^^ since 
?7 ^ 0. 

Inelastic scattering {gg ^ ggg) is implemented in a kinetic approach only recently. 
Although this is a higher order process in perturbative expansion, it turns out to 
affect elliptic flow significantly. See Refs. [79, 80, 81] 

In kinetic theories, momentum exchanges among particles are responsible for 
equilibration. However, forward scattering with very small scattering angle is 
insufficient for the system to equilibrate. So the effective (transport) cross section 
can be defined as (Ttr = / rf^cm sin^ ^cm J)" , where 6'cm is scattering angle in the 
center of mass system between two scattering particles. 

The Boltzmann equation is applied to dilute gases where two particle correlation 
can be ignored. So one should keep in mind the applicability condition of the 
kinetic theory in this case. 
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Through measurement of V2 and its analysis in terms of hydrodynamic/transport 

models, one can extract the transport properties of the matter produced in 
H.I.C. In the next subsection, we discuss hydrodynamic modeling of H.I.C. 

4.3 Ideal hydrodynamic model 

Hydrodynamics introduced in Sec. 3 is a general framework to describe the 
space-time evolution of locally thermalized matter for a given equation of 
state (EoS). This framework has hcni applied to the intermediate stage in 
H.I.C. In this section, we neglect the effects of dissipation and concentrate 
on discussion about ideal hydrodynamic models. The main ingredient in ideal 
hydrodynamic models in H.I.C. is the EoS of hot and dense matter governed 
by QCD. In addition, one also needs to assign initial conditions to the hy- 
drodynamic equations. Hydrodynamics can be applied to a system in which 
local thermalization is maintained. However, in the final state of H.I.C. the 
particles are freely streaming toward the detectors and their mean free path is 
almost infinite. This is obviously beyond the applicability of hydrodynamics. 
Hence we also need a description to decouple the particles from the rest of the 
system. To summarize, the hydrodynamic modeling of H.I.C. needs an EoS, 
initial conditions and a decoupling prescription. Modeling of these ingredients 
in hydrodynamic simulations has been sophisticated for these years and tested 
against a vast body of RHIC data. 

We first look at the EoS in more detail. The EoS is in principle calculated 
from lattice QCD simulations. The realistic results with (almost) physical 
quark masses are obtained recently [53]. However, if one wants to utilize the 
EoS from lattice simulations, one needs to interpret the EoS in term of a 
hadron picture [82] since one calculates momentum distributions of hadrons 
in the final decoupling stage. For this purpose, the lattice EoS is compared 
with the resonance gas model below Tc- If there exists a deviation between 
them, it prevents ones from utilizing the lattice EoS directly in hydrodynamic 
simulations. Instead, in hydrodynamic simulations, the models of EoS depicted 
in Fig. 16 arc conventionally used [45]. The most simple EoS (EOS I) is P = 
e/3 for an ideal gas of relativistic massless particles^^. A more realistic EoS 
(EOS Q) includes the effect of hadron masses and phase transition between 
hadronic matter and the QGP. At low energy density the EoS is described by a 
hadron resonance gas model (EOS H). This particular model includes almost 
all the hadrons in the Particle Data Table [83], while some models include 
only ground states of hadron multiplets or several low mass resonances. At 
high energy density, the EoS can be described by a bag model 

P=\{e-AB). (87) 

This EoS is always obtained in relativistic conformal field theories in which the 
trace of energy-momentum tensor is vanishing T''^ = e — 3P = 0. So the particles 
are not necessarily "free" . 
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Fig. 16. Some typical EoS in hydrodynamic models [45]. 



The bag constant B is tuned to match pressure of the QGP phase to that of 
a hadron resonance gas at critical temperature Tc: Pqgp(T'c) — ^'hadron(?c)- 
As discussed in Sec. 4.1, a hadron gas in H.I.C. is not in chemical equilibrium 
below the chemical freezeout temperature. T'^^ which is closed to Tc, so the 
hadron phase may not be chemically equilibrated in H.I.C. A chemically frozen 
hadron resonance gas can be described by introducing the chemical potential 
for each hadron [84, 85, 86, 87, 88, 89, 90]. The numbers including aU 
decay contributions from higher-lying resonances, Ni = Ni + '^ji bR^ixNjf , 
are conserved during the evolution in co-moving frame of fluid elements. Here 
Ni is the number of the i-th hadronic species in a fluid element and h^^ix 
is the effective branching ratio (a product of branching ratio and degeneracy) 
of a decay process i + X. One can calculate the chemical potential as a 
function of temperature from the following conditions: 



Instead of solving continuity equations for each hadron, the effect of hadron 
number conservation can be embedded in the EoS of resonance gas through 
IJ,i{T) obtained above. 

For a decoupling prescription, the Cooper- Frye formula [91] is almost a 
unique choice to convert the hydrodynamic picture to the particle picture 



s{T,{fi,}) " s{Tc,{^^^} = 0)^ 



(88) 




s 

d f p ■ da{x) 



(89) 



(27r)3 exp [{p ■ u{x) - fiix))/T{x)] ± 1 ' 



(90) 



where E is the energy, / is the phase space distribution, d the degeneracy of 
the particle under consideration {e.g., d = 3 for pions), p is the momentum, da 
is the normal vector to the freezeout surface element, u is the four- velocity, is 
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the chemical potential and T is the decoupling temperature assuming isother- 
mal freezeout hypersurface S. Contribution from resonance decays should be 
taken into account by applying some decay kinematics to the outcome of the 
Cooper-Frye formula. The decoupling temperature T'^'"^ is fixed through si- 
multaneous fitting of pt spectra for various hadrons in the low pT region. In 
the blast wave model, decoupling temperature and radial flow velocity are in- 
dependent parameters to fit pT spectra. On the other hand, there is a negative 
correlation between T'^'^'^ and average radial fiow velocity in the hydrodynamic 
model: the lower decoupling temperature, the larger average radial fiow ve- 
locity. This formula ensures the energy-momentum conservation on freezeout 
hypersurface S as long as the EoS is calculated using the same distribution 
function. If one puts resonances up to the mass of 2 GeV in the resonance 
gas model, one should calculate all the contribution of hadrons in the EoS. 
Otherwise, neglect of the contribution leads to violation of the energy mo- 
mentum conservation^". It should be noted that p ■ da term in Eq. (90) can 
be negative. This means the in-coming particles through S are counted as a 
negative number. Although this seems peculiar, this negative contribution is 
needed for global energy momentum conservation. 

The prescription to calculate the momentum distribution as above is some- 
times called the sudden freezeout model since the mean free path of the par- 
ticles changes from zero (ideal fluid) to infinity (free streaming) within a thin 
layer S. Although this model is too simple, it has been used in hydrodynamic 
calculations for a long time. It is illustrated in Fig. 17(a). Recently one utilizes 



hadronic cascade models to describe the gradual freezeout [92, 93, 94, 95, 96]. 
As will be shown, this hadronic afterburner is mandatory in understanding V2 
data. Phase space distributions for hadrons are initialized below Tc by using 



If the lattice EoS below cannot be described by a resonance gas model, the 
Cooper-Frye formula violates the energy-momentum conservation on S. This is 
the reason why there are only few serious attempts of lattice EoS to hydrodynamic 
simulations. 




(a) Sudden freezeout 



(b) Gradual freezeout 



Fig. 17. Two freezeout pictures in H.I.C. 
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the Cooper-Frye formula. The hadronic cascade models describe the space- 
time evolution of the hadron gas. This model is illustrated in Fig. 17(b). This 
kind of hybrid approaches in which the QGP fluids are followed by hadronic 
cascade models automatically describes both the chemical and thermal freeze- 
out and is much more realistic especially for the late stage. 

Initial conditions in hydrodynamic simulations are so chosen as to repro- 
duce the centrality and rapidity dependences of multiplicity dNch/drj. Initial 
conditions here mean energy density distribution e{x,y,ris) and flow velocity 
u'^{x,y,r]s) at the initial time tq. Again baryon density is neglected since, at 
midrapidity at RHIC, the net baryon density is quite small. The pressure dis- 
tribution can be obtained from the energy density distribution through the 
EoS. Space-time rapidity rjs independent initial energy density distribution 
e{x,y,r]s) = e{x,y) and Bjorken scaling solution Wgj are assumed in (2-1-1)- 
dimensional hydrodynamic simulations. In this case, one discuss the obscrv- 
ables only at midrapidity. At ris = 0, one can parametrize [98] the initial 
entropy density based on the Glauber model 

■^(^' y) = , jo — a«part(2;, y, b) + (l-a)ncoii(x, y; b) (91) 

The soft/hard fraction a is adjusted to reproduce the measured centrality de- 
pendence [99] of the charged hadron multiplicity at midrapidity. By using the 
EoS, one can calculate the initial energy density distribution from Eq. (91). 
For fully three-dimensional initial conditions, see Refs. [84, 94, 97]. A novel 
initial condition is based on the color glass condensate (CGC) picture [100]. 
One can calculate the local energy density of produced gluons within the CGC 
framework [101, 102] and utilize it as an initial condition of hydrodynamic sim- 
ulations. In Fig. 18, an example of the CGC initial energy density distribution 
for a non-central H.I.C. in a full (3-|-l)-dimensional hydrodynamic simulation 
[103] is shown in the transverse plane (left) and in the reaction plane (right). 
In the right side panel the horizontal axis corresponds to the impact param- 
eter direction, and the vertical axis to the space-time rapidity r/s. Figures 19 
and 20 show charged particle multiplicity from hydrodynamic simulations are 
compared with the PHOBOS data [99, 104]. Figure 19 shows dNch/dij as 
a function of the number of participants (iVpart) [99]. This data is fitted by 
using two kinds of initial conditions; from Glauber model calculations and 
from Color Glass Condensate (CGC) model calculations [94]. Both models 
reproduce the centrality dependence of the data. Figure 20 shows the rapidity 
distribution of dNch/drj for each centrality [104]. The fitting of multiplicity is 
the starting point of further analysis based on hydrodynamic simulations. 

In the hydrodynamic models, various combinations of initial conditions, 
EoS and decoupling prescriptions arc available to analyze the experimental 
data in H.I.C. Of course, final results largely depend on modeling of each 
ingredient. So it is quite important to constrain each model and its inherent 
parameters through systematic analyses of the data toward a comprehensive 
understanding of the QGP. 
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Fig. 18. Energy density distribution in a non-central H.I.C. within a CGC initial 
condition in the transverse plane (left panel) and in the reaction plane (right panel). 
The two horizontal thick black lines in the right panel are the Lorentz contracted 
nuclei. The color gradation in the right side of each panel indicates the energy density 
scale in unit of GeV/fm^. 




Fig. 19. Centrality dependence of multiplicity from PHOBOS [104] are fitted by 
hydrodynamic calculations with two different initial conditions [94, 103]. 
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Fig. 20. Pseudorapidity dependence of multiplicity from PHOBOS [104] are fitted 
by hydrodynamic calculations with two different initial conditions [94, 103]. 
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4.4 Application of the ideal hydro dynamic model to H.I.C. 

In this subsection we analyze H.I.C. at RHIC in terms of ideal hydrodynamic 
models discussed in the previous subsection. 

Before we start our main discussion on elliptic flow parameter V2, we 
mention here that the transverse momentum distributions for pions, kaons 
and protons are also important since these reflect dominant transverse flow, 
namely radial flow. Currently, among hydrodynamic models, yields and slopes 
of pt spectra are reproduced in pure hydrodynamic calculations with early 
chemical freezeout or in gradual freezeout approaches. It should be noted 
here that simultaneous reproduction of the yields and the slopes is important. 
Sometimes, one only compares the slope of the pr spectra by scaling the yields 
"by hand" within hydrodynamic approaches. However, chemical composition 
of hadronic matter does affect the transverse expansion [84]. Therefore, it 
does not make any senses if one compares only the slopes by keeping chemical 
equilibrium of hadrons. 

As discussed in Sec. 4.2, V2/S can be interpreted as a response of the sys- 
tem to initial spatial eccentricity. Figure 21 shows V2/e as a function of the 
transverse multiplicity density {1/ S)dNch/dy from AGS to RHIC energies. 
Hydrodynamic results in Fig. 14 are shown symbolically as horizontal lines. 
The experimental data monotonically increase with particle density, while 
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Fig. 21. V2/e as a function of transverse multiplicity density compiled by NA49 
Collaboration [105]. 



ideal hydrodynamic response is almost flat [77]. Ideal hydrodynamics is ex- 
pected to generate the maximum response among the transport models^^ . The 
experimental data reach this limit for the first time at RHIC. Figure 22 shows 

It should be emphasized again that the hydrodynamic results above are obtained 
by a particular combination of modeling, i.e., Glauber type initial conditions, 
EOS Q with chemical equilibrium in the hadron phase and sudden freezeout at 
fixed decoupling temperature. 
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the differential elliptic flow, f 2 as a function of transverse momentum for pi- 
ons, kaons, protons and lambdas. A mass ordering pattern is seen in V2 data, 
which was predicted by ideal hydrodynamic calculations [106] The pseudo- 



012I- I < -i.<j 

(a)200 GeVAu + Au y'.- yy 

(mininmm bias) y' 




0.2 0.1 O.e 0.8 1.0 1.2 1.4 1.6 

Transverse momentum p-, (GeV/c) 



Fig. 22. Differential V2 for pions, kaons, protons and lambdas [64]. 

rapidity dependence of observed by PHOBOS [107] has a triangular shape 
as is seen in Fig. 23. In the pure ideal hydrodynamic result, hydrodynamic 
equations are initialized by the Glauber model and are solved all the way 
down to T^'^'^ = 100 MeV. The pure hydrodynamic model gives a comparable 
result with the data only at midrapidity. However at forward and backward 
rapidities, it overshoots the data significantly. If we replace the hadron fluid 
with a hadron gas utilizing a hadron cascade, V2 is significantly reduced in the 
forward and backward region. In this hybrid model the hadrons have a finite 
mean free path, which results in an effective shear viscosity in the hadron 
phase. So dissipative hadronic "corona" effects turn out to be important in 
understanding the V2 data. The model also reproduces a mass ordering pattern 
of vi for identified hadrons as a function of near midrapidity in Fig. 24. 
Figure 25 shows the centrality dependence of ■ The solid line is the result 
from ideal hydrodynamic calculations, while the dotted line from the hybrid 
model. It is clear that for peripheral collisions, where the multiplicity is small, 
the hadronic viscosity plays an important role. One may notice that the re- 
sult from the hybrid model systematically and slightly smaller than the data. 
However, there could exist the effect of initial eccentricity fluctuations which 
is absent in this hydrodynamic calculations. The deviation between the results 
and the data can be interpreted quantitatively by this effect. Figure 26 shows 
V2{pt) for pions, kaons and protons in 10-50% centrality at 77 = (left), t] = 1 



There is a caveat to interpret the agreement since this particular hydrodynamic 
calculation does not reproduce particle ratios due to a lack of early chemical 
freezeout. The importance of hadronic viscosity and chemical freezeout in hydro- 
dynamic calculations is recognized [108] after the announcement of the discovery 
of perfect fluid QGP [2]. 
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Fig. 23. Pseudorapidity dependence of V2. PHOBOS data [107] compared to differ- 
ent model calculations [94]. 
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Fig. 24. Differential V2- STAR data [64] compared to model calculations [95]. 



0.14p 
0.12: 
0.1: 
0.08- 

■1 

0.06 
0.04 
0.02 






T,,„=100MeV 






hydro+cascade 




• PHOBOS{hit) 






• PHOBOS(track) 




,.-{|-i.t.!4>i^ 












. 1 .... 1 .... 1 .... 1 .... 1 .... 1 . 


I 



50 100 150 200 250 300 350 400 



Fig. 25. V2 as a function of centrality. PHOBOS data [107] compared to different 
model calculations [94]. 
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(middle) and 77 = 3 (right) observed by BRAHMS [109]. Also here the hybrid 
model reproduces the slope of these differential elhptic flow parameters. 




We would like to point out here that the mass ordering, clearly visible 
in Fig. 24, is there in the final result. If one would look at the result just 
after the QGP phase transition, the difference between the pions and protons 
would be quite small. So it turns out that the splitting patterns are caused 
by hadronic rescattering. This is illustrated in Fig. 27. One can conclude that 
the large magnitude of the integrated V2 and the strong mass ordering of the 
differential V2{pt) observed at RHIC result from a subtle interplay between 
perfect fluid dynamics of the early QGP stage and dissipative dynamics of 
the late hadronic stage: The large magnitude of V2 is due to the large over- 
all momentum anisotropy, generated predominantly in the early QGP stage, 
whereas the strong mass splitting behavior at low px reflects the redistribu- 
tion of this momentum anisotropy among the different hadron species, driven 
by the continuing radial acceleration and cooling of the matter during the 
hadronic rescattering phase. 

We have seen so far that the hydrodynamic model which includes Glauber 
type initial conditions followed by a perfect fluid QGP and a dissipative 
hadronic gas evolution is the most successful combination for describing the 
RHIC data. We now go to the discussion on the initialization dependence of 
V2- Two types of initial conditions, namely the Glauber type initial conditions 
and the CGG initial conditions, are discussed in the previous subsection. V2 
as a function of centrality is shown again for these two initial conditions in 
Fig. 28. In the case of the Glauber initial conditions we can conclude early 
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Fig. 27. Differential V2 with and without hadronic rescattering [95]. 

thermalization and the discovery for the perfect fluid QGP. In the case of 
the CGC initial conditions, we cannot, however, claim the discovery since 
the model initialized by CGC overshoots the data in almost the whole range. 
Since the hydrodynamic model calculations depend on the initial conditions, 
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Fig. 28. V2 as a function of centrality. PHOBOS data [107] are compared to hydro- 
dynamic results with two different sets of initial conditions [94]. 

it is very important to understand them before making final conclusions. In 
the case of CGC initial conditions viscosity might be needed even in the QGP 
phase to get the model down to the data points. The effect of viscosity could 
therefore be quite important. The high V2 values from the CGC initial con- 
ditions are traced back to the initial eccentricity. In Fig. 29(a) the energy 
density distribution in the impact parameter direction is plotted for different 
conditions. If the energy density profile has a sharp edge (no diffuseness), an 
integral in Eq. (78) is relatively weighted in the edge region and, consequently, 
eccentricity becomes maximum at a given impact parameter. If one compare 
the energy density profile of the CGC with the one of the Glauber model, one 
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see the CGC profile has a sharper edge than the Glauber model does. The re- 
sultant eccentricity as a function of impact parameter is shown in Fig. 29(b). 
Eccentricity from the CGC is about 20-30% larger than that from the Glauber 
model. This is the reason why hydro + hadronic cascade approach which even 
includes hadronic viscosity overshoots the V2 data 




X (tm) b (fm) 

(a) Energy density distribu- (b) Initial eccentricity 

tion 

Fig. 29. The difference eccentricity between Glauber and CGC initial conditions. 



4.5 Summary 

Hydrodynamics is a framework to describe the space-time evolution of matter 
under local equilibrium. It is applied to the intermediate stage in H.I.C. to 
extract the transport properties of the QGP from RHIC data. Hydrodynamic 
modeling includes initial conditions, EoS and decoupling prescriptions. Final 
results certainly depend on combination of each modeling. So much attention 
should be paid to these ingredients before drawing robust conclusions from 
hydrodynamic analyses. Elliptic flow has played a mayor role in understanding 
the transport properties of the QGP. Glauber initial conditions, ideal hydro- 
dynamics in the QGP phase and dissipative gas for the hadron phase are three 
pillars for agreement between the model and elliptic flow data. Whereas, if 
CGC initial conditions are employed, the initial eccentricity gets increased by 
20-30%. If the nature chooses this kind of initial condition, viscosity might be 
needed even in the QGP phase. 
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